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Abstract 

We calculate the contribution of the fluctuations with the thermal origin to 
the inflationary nonGaussianity. We find that even a small component of radi- 
ation can lead to a large nonGaussianity. We show that this thermal nonGaus- 
sianity always has positive /nl- We illustrate our result in the chain inflation 
model and the very weakly dissipative warm inflation model. We show that 
Inl ~ is general in such models. If we allow modified equation of state, 
or some decoupling effects, the large thermal nonGaussianity of order /nl > 5 
or even /nl ~ 100 can be produced. We also show that the power spectrum 
of chain inflation should have a thermal origin. In the Appendix A, we made 
a clarification on the different conventions used in the literature related to the 
calculation of /nl- 
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1 Introduction 



Inflation has been remarkably successful in solving the problems in the standard hot 
big bang cosmology [H [21 [31 H] . Furthermore inflation shows us that the fluctuations 
of quantum origin were generated and frozen to seed the wrinkles in the cosmic 
microwave background (CMB) [51 [6] and today's large scale structure [TJ, [HI El Ell 
[TT] . Its prediction of a scale invariant spectrum which has been confirmed in the 
experiments in the past decade is remarkable and has been taken to be a great success 
of the theory. 

Since the idea of inflation was proposed, there have been a large number of infla- 
tion models. It has become one of the key problems in cosmology to extract more 
information from experiments in order to distinguish these inflation models. The 
key quantities from experiments include the power spectrum of scalar and tensor 
perturbations, the scalar spectral index and its running, and the nonGaussianity. 

The amount of nonGaussianity is often estimated using the quantity /nl, which 
can be written a£ 

c = c 9 + jj/NL (c 9 2 - <0) , (1) 

where the subscript g denotes the Gaussian part of £. 

It has been shown that in the simplest single field slow roll inflation models, 
the nonGaussianity estimator / NL ~ 0(e,r)') [131 [33], where e and rj are the slow 
roll parameters. Such a small nonGaussianity is not only much smaller than the 
current observational bound |/nl| < 100, but also well below the sensitivity of the 
Planck satellite, /nl ~ 5. However, the recent study of the inflation models with 
significantly nonlinear dynamics shows that the nonGaussianity in them could be 
large. For example, in the DBI [15], K-inflation [16], and ghost inflation [!17] models, 
/nl can reach 0(1) or larger than 0(1) in some parameter regions. Such a large 
nonGaussianity is hoped to be observed in the future experiments and shed light on 
the physics behind these models. 

Recently, Yadav and Wandelt has claimed that from the WMAP 3-year data, /nl 



1 Note that for the definition of /nl, there is a sign difference between the notation of the WMAP 
group [123 an d Maldacena's calculation [T3] . We use the same notation as that of the WMAP group 
here. The difference is discussed in detail in the Appendix A. 
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is detected at above 99.5% confidence level [18]. They show that at 95% confidence 
level, the local shape /nl is in the region 

26.91 < /nl < 146.71 . (2) 

If this result is confirmed by the WMAP 5-year data and Planck, a great number of 
inflation models (without extra mechanisms) will be ruled out, and there is also hope 
to measure the shape of /nl, and the tri-spectrum tnl at Planck. 

In this paper, we propose another mechanism to produce potentially large non- 
Gaussianity. Instead of producing nonGaussianity from the nonlinear evolution of the 
inflaton, in our mechanism, the large nonGaussianity stems from the correlation in 
the initial conditions. We will show that if the initial condition of the perturbations 
is prepared in part by thermal fluctuations, there can be strong 3-point correlation, 
inducing large nonGaussianity. 

In many inflation models, the radiation component only takes a very small part 
in the energy density. But since the nonGaussianity from the coherent motion of 
inflaton is highly suppressed, the thermal nonGaussianity can play a significant part 
in /nl, and in some parameter regions it provides the dominant contribution. This 
may open a window for us to study the thermal fluctuations in the models. 

Indeed, the thermal effects are significantly important in some inflation models. 
One example is the chain inflation. Based on the the rapid tunnelling mechanism for 
the meta-stable vacua in the string landscapeJTH [20], Freese, Spolyar and Liu [21] 
proposed the so-called chain inflation model in which the meta-stable vacua during 
inflation tunnel very rapidly. The density perturbation in chain inflation is calculated 
by Feldstein and Tweedie in [22] and a simplified version of the chain inflation was 
proposed by Huang in |23j. 

In the chain inflation models, the average life time for a meta-stable vacuum is 
much smaller than the Hubble time, so that the vacuum decay via bubble nucleation 
takes place very rapidly, and there can be many bubbles nucleated within one infla- 
tionary horizon. These bubbles eventually collide and the energy stored in the bubble 
wall decays into the radiation. This is very different from the slow roll inflation models 
in which the decrease of the inflaton energy density is wasted by the cosmic fraction, 
with very little radiation being left. 
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Another example with large thermal effect is the warm inflation by Berera and 
Fang [24]. In the warm inflation model, a fraction of inflaton energy decays into 
radiation continuously during inflation. The decay from inflaton to radiation can be 
achieved by a interaction term in the inflaton's Lagrangian. It is shown in warm 
inflation that due to the continuously creation of radiation, the temperature during 
inflation can be nearly constant [23] , so it provides a playground for investigating the 
thermal effects. Previously, the nonGaussianity of the warm inflation model is studied 
in [25]. But in [25 j . the authors considered only the nonGaussianity of the inflaton 
field with Gaussian noise source, and the nonGaussianity of the thermal fluctuation 
has not been investigated. In this paper, we only consider the warm inflation in the 
very weakly dissipative regime, in which the existence of the thermal bath would not 
spoil the quantum vacuum. 

In such inflation models with thermal radiation, it can be shown that the non- 
Gaussianity estimator /nl is no longer suppressed by the slow roll parameters. Even 
when the radiation component is so tiny that it does not qualitatively change the 
inflationary background, considerable nonGaussianity /nl ~ 0(1) can be produced. 

Furthermore, we suppose in some cases, a new scale related to the acoustic horizon, 
or some decoupling scales enters the calculation. In this case, very large nonGaus- 
sianity of the order /nl > 5 or even /nl ~ 100 can be produced without fine-tuning. 

This paper is organized as follows. In Section 2, we develop the general method 
to calculate the nonGaussianity of the thermal origin. We calculate the 2-point and 
3-point correlation functions of thermal fluctuations. Based on these, we derive the 
power spectrum and the nonGaussianity estimator /nl- In Section 3, we calculate 
the amount of nonGaussianity explicitly in the chain inflation model and the thermal 
inflation model. We conclude in Section 4. 

2 The thermal correlation functions and the non- 
Gaussianity 

In this section, we calculate the correlation functions, the power spectrum and the 
nonGaussianity of thermal fluctuations. We also give a simple estimate of the non- 
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Gaussianity by calculating the backreaction. 

We suppose the energy density takes the form 



p EE po + Pr = PO + AT m , (3) 

where po is the energy density without thermal origin, for example, the effective 
vacuum energy provided by the inflaton potential. And p r is the energy density for 
the radiation, and A is a constant with dimension [mass] 4 ~ m . Note that m = 4 for 
usual radiation. While for generality, phenologically, we still keep m here. 

The correlation functions in thermal equilibrium can be calculated from the par- 
tition function of the system 

Z = XV^, (4) 

r 

where (3 — T -1 . 

Let U ee pV represents the total energy inside a volume V. Then the average 
energy of the system is given by 

(U) = -^f , (5) 

The 2-point correlation function for the fluctuations Sp = p — (p) is given by 

2 _ (5U 2 ) 1 d 2 \ogZ 1 _ mAT m+1 

( P } ~ V 2 ~ V 2 dp 2 ~ V 2 d(3 V ' U 

where in the final equality we have neglected "()" because the difference is next to 
the leading order. 

Similarly, the 3-point correlation function can be expressed as 

3 _ (5U 3 ) _ 1 dHogZ 1 d 2 (U) _ m(m + l)AT m+2 

\ 5 p i - -yr- ~ d p3 ~V^~dp 2 ~~ V 2 ' () 

Now let us apply the above calculation to inflation. First, we calculate the equation 
of state w r for general radiation p r = AT m . For simplicity, we only consider the case 
that w r is a positive constant. To do this, we temporarily consider radiation without 
source. In the expanding background, consider a comoving volume V c which is in the 
thermal equilibrium. The conserved radiation entropy within this volume is given by 

S = P -^V C . (8) 
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The radiation energy density and pressure changes with respect to the scale factor a 
as 

Pr ~ Pr ~ a- 3(M, ' +1) . (9) 

Combining (jSJ) and ([9]), the temperature scales as 

T~(p r +p r )V c ~a- 3w r , (10) 
so the relation between p r and T can be written as 

p r ~ T «"r . (11) 

So we have the relation between the sound speed, the equation of state, and the 
parameter m defined in (J3J) as 

C 2 s =W r = . (12) 

m — 1 

Another important issue is to determine the appropriate size of the thermal system 
L. By determining L, we mean that at length scales smaller than L, the fluctuation 
of the system can be calculated using the thermal dynamics described above, and at 
scales greater than L, the fluctuation is governed by the cosmological perturbation 
theory. Note that a typical photon in the thermal system has wavelength T _1 , so 
there is a lower bound L > T _1 on L. Otherwise, the system is too small to be 
treated as a thermal system, and the above calculation no longer holds. Also, there 
should be no thermal correlation outside the acoustic horizon c s i7 _1 , so the constraint 
on L is T- 1 < L < c s H~\ [» 

We will argue in the discussion that from some decoupling mechanism, the ex- 
plicit value of L may depend on the detailed properties of the thermal system and 
the dynamics of inflation. While in the remainder of the paper, we will hold L as 
a parameter (sometimes called the "thermal horizon") during the calculation, and 
discuss the most modest limit L = c*H~ l when we come to final results. 



2 In the literature, the length scale T -1 is used as L to calculate the thermal fluctuations by some 
authors |26j . In our calculation of nonGaussianity, if we choose L = T^ 1 , the result turns out to be 
much more dramatic: we will get very large nonGaussianity for a much wider class of inflationary 
models. While in our paper, we do not choose to use L = T _1 , and only take T~ 1 as an lower bound 
of L here. 
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Figure 1: This figure illustrates how the initial condition of perturbation is prepared by 
the thermal fluctuations. The black cycle represents the thermal horizon. A fluctuation 
Sp of the system can exit the thermal horizon L (shown as a shell in the figure) during 
inflation. This shell outside the thermal horizon can not return to thermal equilibrium with 
respect to the original volume. It provides the initial condition of inflationary fluctuations. 
Fluctuations are created shell by shell. 

The fluctuation Sp can be thought of as an average over all the local fluctuations 
of the thermal system 5p = y J d 3 xSp(x). Performing the fourier transformation, 
and linking the zero mode of k to the horizon exit mode k = a/L, as illustrated in 
Fig. [TJ Then the out-of-thermal-equilibrium initial condition for Sp k takes the form 



As a check, the equation ffl3|) can also be obtained using the window functions. 

Since L is smaller than the inflationary horizon, the relation between the energy 
density perturbation and the scalar type metric perturbation at the boundary of V 
can be calculated using the Poisson equation [] 



where G is the Newton constant, and is the fourier mode of the Newtonian 
gauge metric perturbation defined as ds 2 = a 2 (—(1 — 2&)dr] 2 + (1 + 2$)dx 2 ), and is 
calculated at k = a/L. Note that $ here is not the Newtonian potential <3>at, but 

3 Note that to be exact, the 00 component of the linearized Einstein is — Vp h $ + 3if $ + 3H 2 <& = 
AitGSp. The sound speed c s do not enter this equation. So as long as c s is not too large, even when 
L = CgH^ 1 , the equation (Tl4|) is a very good approximation. 



5p k = k 2<5p . 



(13) 



$ fcL = 4nG5p k L 2 , 



(14) 
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rather $ = — $at. We follow the WMAP convention to use $ as the perturbation 
variable. Further discussion on the conventions can be found in the Appendix A. 

We work in the Newtonian gauge only for simplicity. Since L < iif -1 , the different 
choice of gauge is not important inside the thermal horizon (see, for example, |27j). 
After the inflationary modes leave the thermal horizon, we use $ to describe the 
mode. It is well known that $ can be made gauge invariant when considering the 
more general gauges. So our final result will be independent of gauge choice. 

The equation ( TT4l) provides a thermal initial condition of After that, the 
evolution of is governed by the cosmological perturbation theory. In the lowest 
order slow roll approximation, ^ evolves as 

$k ~ V^H^-kr) , (15) 

where r is the comoving time and H^ 2 is the first kind Hankel function. It can be 
shown that $fc oscillates inside the inflationary horizon with nearly constant ampli- 
tude, so \<&k\k=aH — \&kh\- After horizon crossing, the amplitude of is frozen so 
that the change of during a few e- folds is negligible. Using ([6]), ffT3l) and ffT4l) . the 
2-point correlation for at k = aH (and also a few e-folds outside the inflationary 
horizon) is expressed as 

($2) = (47rG) 2 mp r TLk- 3 . (16) 

Note that here k = |k|, and in (1T61) . we are actually calculating the correlation between 
the mode k and — k. Similarly, in the three point correlation function, the quantity 
we calculate corresponds to the equilateral triangle, satisfying k x + k 2 + k 3 = and 
|ki| = | lt2 1 = [IC3 1 = k. The phase of in the correlation functions cancels due to 
the momentum conservation. 

The power spectrum of from the thermal origin can be written as 

P^^(<S>l} = 8G 2 m Pr TL. (17) 

Similarly, the 3-point function of $fe at k = aH is expressed as 

($|) = (4vrG) 3 m(m + l)p r T 2 k~ 9/2 . (18) 

The nonGaussianity can be calculated for (k, relating to by ( ~ $/t/e at a few 
e-folds outside the inflationary horizon, where e = —H/H 2 . 
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Note that this 3-point correlation function (Cf ) is always positive. The positivity 
is transparent when we choose L as the acoustic horizon c s i/ _1 . In this case, C freezes 
outside L and do not have chance to change sign. While in the case that L < c s H~ l , 
although ( has oscillating solution inside the acoustic horizon, but (Cf) still can not 
change sign. This is because after the initial condition is prepared, the evolution of 
(Cf ) is governed by the interaction Hamiltonian of C as 



We are assuming the slow roll inflationary scenario, and do not employ other mecha- 
nisms to have large interaction H int . So once the positive nonGaussian initial condi- 
tion is produced, it will keep positive until it is observed in the CMB. 
Finally the nonGaussianity estimator /nl takes the form 



Note that here we have assumed that the origin of perturbation is completely ther- 
mal. A combination of the thermal and quantum origin of the power spectrum and 
nonGaussianity will be discussed at the end of this section. 

From ( l20~j) . we see that /nl oc L~ 2 . So the smaller L is, the larger the nonGaus- 
sianity can be. Note that for modified equation of states, if \m\ <C 1, the m in the 
denominator also enhances the nonGaussianity. 

The nonGaussianity (1201 could also be estimated without calculating the 3-point 
correlation function explicitly. The idea is to calculate the back-reaction. A fluctua- 
tion mode which crosses the thermal horizon earlier can change the background for a 
later fluctuation mode, so leads to nonGaussian correlation between these two modes. 

From the 2-point correlation function, or from the standard result in thermody- 
namics that 5T/T ~ a/1/CV, we have for the second mode 



As the first mode has crossed the thermal horizon by the time the second mode 
crosses the thermal horizon, the first mode leads to a modification of the background 
of the thermal system for the second mode. This modification of the background 




(19) 




18 (C|) (C|) ^Gmp r I? ' 



(20) 



5 2 p ~ y/mAT m+1 /V . 



(21) 
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represents the correlation of the two modes, so is the nonGaussian contribution. At 
the thermal horizon, this nonGaussian contribution takes the form 



'" + i^T.-VW~<^. (22) 



2 v ' 2V 

As discussed earlier, when the modes reaches the inflationary horizon k = aH, we 
have 

C - C, ~ 4nG5pL 2 /e . (23) 

So finally the nonGaussianity /nl is estimated as 

5(C-C 3 ) 5e(m+l) 
/NL ~ —3(f- ~ 24^^ ' (24) 

This differs from ( 120]) only by a factor of 3, and can be considered as in good agreement 
for a rough estimate. This back-reaction estimate provides a check for the 3-point 
function calculated above, and also explains why the nonGaussianity can be so large: 
the thermal horizon is smaller than the inflationary horizon, so a back-reaction cal- 
culated at the thermal horizon is larger than the one calculated at the inflationary 
horizon. This large back-reaction leads to a large nonGaussianity. Note that although 
this estimate can not give the precise shape of /nl, the limit we take is similar to the 
squeezed limit, which leads to a local shape nonGaussianity. 

Generally, there can also be perturbations from the vacuum fluctuations of the 
coherent rolling inflaton field. Let us denote this perturbation by $™ c . Since the 
vacuum fluctuation and the thermal fluctuation are of the different origin, they do 
not have correlations between each other. So in the 2-point and 3-point functions, 
the cross terms such as ($fc$^ ac ) vanishes. So for the total power spectrum and the 
nonGaussianity, 

ptot = pvac + p f tot = O ^ _^ / i- ^ k f 

18 {( Cr 2) + <0} 2 ' 

where Q ac and P,| ac are the comoving curvature perturbation and the power spectrum 
calculated from the vacuum fluctuation of the inflaton field. When |$™ c | <C \&k\, (ESI) 
returns to (ITT)) and (120)) . and when |$™ c | ^> |$fc|, (125)) returns to the power spectrum 
and the nonGaussianity with zero temperature. 
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3 Examples of inflation models with large thermal 
nonGaussianity 



In the previous section, we have given the general formalism to calculate the thermal 
perturbations and nonGaussianity. In this section, we apply the formalism to the 
chain inflation and the warm inflation. 

Although the radiation energy density can be inflated away very easily, as dis- 
cussed in the introduction, there are mechanisms to continuously produce radiation 
so that the radiation energy density keeps nearly constant. Such mechanisms include 
the interaction of the radiation with the inflaton, and the bubble collision during the 
chain inflation, which we shall show explicitly. 

In the chain inflation model, the vacua tunnel rapidly and the time evolution of 
the vacuum energy density can be approximated by 

p (t) = p (0) - at , (26) 

where a denotes the averaged decay rate of the vacuum energy, (a = - in the notation 
of [23]). Suppose that the decreasing energy converts completely into the radiations 
through bubble collision. Taking into consideration of the red shift of the radiation 
during inflation, the radiation energy density satisfies 

dp r (t) = adt — 3H(1 + w r )p r dt . (27) 

By taking the stationary limit t 3> [(1 + w r )H]~ l , we have 

a 2ep _ H 4nG a 

Pr = 3{1 + w r )H = 3(1 + w r ) ' C = ~lP = ' ( ' 

The relation fl28l) could also be obtained from assuming that the radiation density is 
produced within one Hubble time, then from fl26l) and taking t ~ H~ x , we get directly 
that p r is of the order a/H . 

Note that p r is a slow roll quantity during inflation. As p r ~ T m , when m is not 
too small, T also changes very slowly during inflation. This verifies the assumption 
that the radiation density and the temperature during inflation are almost constants. 

Several mechanisms have been proposed to calculate the fluctuations in the chain 
inflation model. In [22], the authors showed that the perturbations can come from 
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different pathes along which the meta-stable vacuum tunnel. In [23], the density 
perturbation is calculated by applying the standard formulism to the effective scalar 
field characterizing the tunnelling effect. But up to now, the exact mechanism how 
the perturbation in chain inflation is generated is still not clear. 

In this section, we propose a new mechanism to produce the density perturbation 
of the chain inflation. We claim that the perturbation can have a thermal origin. It 
is because after the bubble collision, the energy contained in the bubble wall becomes 
the radiation. As to be shown later in this section, the radiation density is of order 
p r ~ ep, whose thermal fluctuation can exit the horizon and produce a scale invariant 
power spectrum. So in a realistic calculation, the fluctuations discussed in [22], [23] 
and the thermal fluctuation should be taken into account at the same time. 

Using ( [171) . the power spectrum takes the form 

16G 2 mepTL _ 16G 2 mpTL 

P *~ 3(1 + ' C " 3(l + w r )e • ( 9) 

Assuming m and w r are exactly constants, the spectral index takes the form 



n e — 1 = 



dlnP, 



din k 



If L saturates its lower bound L ~ T -1 , and considering the usual type of radiation 
m = 4, w r = 1/3, then we recover the power spectrum and the spectral index in the 
simplified chain inflation model |23j . 

If the thermal perturbation is dominate over other perturbation sources during 
chain inflation, then the nonGaussianity can be read off from (120!) that 

5(m + l)(l + w r ) 
lSm{LH) z 

Note that for the usual type of radiation m = 4, w r = 1/3, L should satisfy L < 
(3i?) _1 , so /nl is always larger than 0(1). This is very different from the ordinary 
inflation model in which /nl ~ ^( e )- 

To go one step further, let us consider the modified radiation m/4. And make 
a modest estimate that L = c s H~ l . In this case, 

5(m + l)(l + w r ) 5(m+ 1) 

Jnl — tt. 2 — To \ 6l > 

ISmcj 18 

As a phenological model, if m is large, one can have large /nl- 
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If the thermal perturbation is not the dominate source in the power spectrum, 
then we need to compare the thermal and other contributions to estimate the non- 
Gaussianity. Let the power spectrum from the other origin be P$ ac , then using (1201) 
and ( l25l) . where we have neglected the nonGaussianity produced by the sources other 
than thermal, we can express /nl as 

_ 5(m + l)(l + w r ) 

JNL — — • 1<J<JJ 

18m{LH) 2 {l + 

So it is clear that although the nonGaussianity is suppressed by the ratio of the 
spectrums, there is no longer the e suppression. Moreover, if L <C H^ 1 , then the 
nonGaussianity can be enhanced by a great amount. 

In the case of warm inflation, the radiation is continuously produced during slow 
roll inflation. This process can be modeled by adding a interacting term between 
inflaton and radiation component in the Lagrangian. In the slow roll regime, the 
equation of motion for inflaton <p is 

3#<p + F v (p + d v V((p) = , (34) 

where Y v is the decay rate for the inflaton to radiation process. We assume T v is a 
constant (or at least a slow roll quantity) here. The equation for the radiation energy 
density takes the form 

p r + 3H(1 + w r )p r = T^if 2 . (35) 
A solution for these equations is given by 

3H{1 + w r ) 

And in [24] . it is shown that this solution is an attractor solution, and independent 
of the initial conditions for the thermal component. 

When p r < p , the universe accelerates. To give a nearly scale invariant spectrum, 
we require the universe undergoes a quasi-dS expansion, so the bound for the radiation 
energy should be p r < epo- This corresponds to the case that T^, < H. When this 
bound saturates, the power spectrum and the nonGaussianity coincides with the chain 
inflation case. 

It can be checked that when the constraint p r < epo is satisfied, the e-folding 
number and the slow roll condition is qualitatively the same as the p r = case, and 
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the inflaton vacuum is not thermalized because the interaction rate < H. So the 
inflationary background and the amplitude of inflaton fluctuation do not change. In 
[21], the authors also considered the case p r > ep , but we will not investigate this 
case in detail here, and assume p r < epo from now on. 

Similarly to the chain inflation case, the thermal power spectrum of warm inflation 
takes the form 

P$ = 8G 2 mp r TL , P c = 8G 2 mp r TL/e 2 . (37) 
And the nonGaussianity of the warm inflation is 

5(m + l)ep 

JNL — , {do) 

27mp r (L#) 2 (l + ^) 2 

which behaves like the chain inflation model with sub-dominate thermal fluctuation. 
Note that p r < epo, so there is no 0(e) suppression. The enhancement due to (LH)~ 2 
still exists. 

To see how /nl behaves at small p r limit, taking into consideration the observa- 

c 



tional constraint of the total power spectrum P} ot ~ 2.5 x 10 9 , /nl can be written 



/nl ^™ • (39) 



as 

40m(m + l)G 3 p r T 2 

DtC 

So when taking the p r — > limit, we get /nl — >• 0. It is reasonable because the /nl 
we are considering comes from thermal origin. 

It seems surprising at first sight that in the limit that e is very small, the radiation 
component bounded by p r < ep is tiny, but the thermal nonGaussianity can still be 
quite large. The reason for this is that given the amplitude of the observed CMB power 
spectrum, when e is very small, the primordial density fluctuation needed to generate 
the CMB power spectrum also becomes small. So a tiny part of radiation becomes 
comparable with the inflaton vacuum fluctuation in the function of generating the 
fluctuations. 

4 Conclusion and discussion 

As a conclusion, in this paper, we have calculated the nonGaussianity from thermal 
effects during inflation. We calculated the 2-point and 3-point thermal correlation 
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functions, and using these correlation functions to calculate the scalar power spectrum 
Pj> and the nonGaussianity estimator / NL . We also used an independent method to 
check the value of /nl- 

We have applied our treatments in the chain inflation. We find that the density 
perturbation in the chain inflation may come from the thermal fluctuations. This 
provides a candidate for the origin of power spectrum of the chain inflation. 

We also calculated the nonGaussianity of the chain inflation model. We found if 
the thermal perturbation is the main source of chain inflation, then the nonGaussian- 
ity /nl of chain inflation is greater than 0(1). Taking into consideration the modified 
sound speed, the nonGaussianity can become much larger. 

If the thermal perturbation is sub-dominate during chain inflation, then / NL is 
suppressed by Pj/(P$ + P| ac ) 2 . But still, there is no 0(e) suppression, and the term 
(LH)~ 2 can provide a large nonGaussianity. 

As another application, we studied the nonGaussianity in the warm inflation 
model. The result of the warm inflation model is similar to the case of chain in- 
flation with sub-dominate thermal component. 

We only studied the p r < ep case in the warm inflation scenario. It is shown that 
large nonGaussianity can already show up in this case. We have not considered in 
this paper the complementary case p r > ep . In this latter case, the thermalization 
of the inflaton vacuum dominates the power spectrum. The thermal part of the 
nonGaussianity /NL Crmal is suppressed in this case. But the thermal nonGaussianity 
of the inflaton field should be taken into consideration. 

In the warm inflation case, it is clear that the vacuum fluctuation and the thermal 
fluctuation are two different sources of inflation fluctuations. So it may lead to large 
isocurvature perturbation. This isocurvature perturbation issue is not discussed in 
detail in this paper. 

In this paper, we calculated the equilateral shape nonGaussianity. And in the back 
reaction estimate, the shape is something like local shape. Since the nonGaussianity 
from thermal fluctuations can be large, and is very hopeful to be observed in the near 
future, it is also important to calculate the more general correlation functions with 
arbitrary k, and obtain the shapes of the nonGaussianity. 

Another important issue is to determine the parameter L given a inflation model, 



15 



which requires more details on the dynamics of the system. In this paper, we mainly 
discussed the upper bound of L, which is governed by the sound speed c s of the 
radiation component. But we note that L may be much smaller than the acoustic 
horizon. One possible mechanism generating smaller L is the decoupling of the fourier 
mode of the thermal fluctuations. When the universe expands, the interaction rate T 
for the thermal fluctuation fourier mode may decrease, so it decouples before reaching 
the acoustic horizon. We wish we will address this issue in the near future. 

The generalization of our calculation to other inflation models with radiation is 
straightforward. For example, our calculation can also be applied to the thermal 
inflation model [28J, or the thermal version of the noncommutative inflation model 




The similar analysis can also be performed in the string gas model [29] , where the 

power spectrum also has a thermal origin. The calculation of the nonGaussianity of 
the string gas model will be represented in a separate publication [30] . 
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Appendix A: Clarification on Conventions 

In this appendix, we clarify the convention we use. To write an appendix to clarify the 
convention is necessary, because a confusion in the convention (especially the sign) 
can lead to an extra minus sign in /nl, and lead to completely opposite predictions. 
This is very different from the calculation of the power spectrum, where a confusion 
of the sign convention usually leads to the same result. 
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In this paper, we use the WMAP convention. In this convention, the metric 
perturbation $ can be written as (in the Newtonian gauge) 

ds 2 = a 2 (-(1 - 2<$>))drf + (1 + 2$)rfx 2 ) . (40) 

So <3> is not the Newtonian potential <3>at, but rather $ = — $jy. This is the same 
as the convention used in [31], which is also the same as the convention used in the 
WMAP group. One can refer to [31J to find the complete definitions. This is of 
the different sign from the convention $ used in [27]. (In [27], they use </> to denote 
the Newtonian potential, and the gauge invariant quantity $ = in the Newtonian 
gauge.) 

For the quantity £, there are also different conventions in the literature. In this 
paper, following the convention of [51] . ( can be written as 

C = $ - ^-<V , (41) 

where (p and 5ip are the background value and the perturbations for the inflaton field 
respectively. And in [27], the ( parameter they use is of the different sign. 

One simple way to check the sign is to relate it to the quantities which have clear 
physical meaning. There are at least two such quantities: the energy density 5p and 
the CMB temperature fluctuation AT/T. In our convention, the Poisson equation 
takes the form 

V 2 

_ $ = 4vrG<5p . (42) 
or 

And the CMB temperature fluctuation can be written as 

7-5— 

Although not related to this paper, we also would like to remind the reader two 
more differences in the conventions, which may be used in the calculation of /nl- One 
is that in [13], Maldacena uses the same convention of ( as that of the WMAP group, 
but the equation in the footnote 16, ( = — 1$, does not follow the WMAP convention. 
So the /nl defined in [13] is of the different sign from the WMAP convention. The 
other is that for the so called "comoving curvature perturbation" TZ. The TZ used in 
[3T] (in the comoving gauge) is of the different sign from the TZ used in [32] outside 
the horizon. 
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